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The use of the ellipsometer for the measurement of the thickness and refractive index 
of very thin films is reviewed. The Poincare sphere? representation of the state of polariza- 
tion of light is developed and used to describe the reflection process. Details of the operation 
of the ellipsometer are examined critically. A computational method is presented by which 
the thickness of a film of known refractive index on a reflecting substrate of known optical 
constants may be calculated directly from the ellipsometer readings. A method for comput- 
ing both the refractive index and thickness of an unknown film is also developed. These 
methods have been applied to the determination of the thickness of an adsorbed water layer 
on chromium ferrotype plates and on gold surfaces. In the former case the thickness was 
23 to 27 A, and in the latter was 2 to 5 A. The measurement of the thickness and refractive 
index of barium fluoride films evaporated on chromium ferrotype surfaces is used as an 
illustration of the simultaneous determination of these two quantities. 



1. Introduction 

Ellipsometry is a convenient and accurate tech- 
nique for the measurement of thicknesses and 
refractive indexes of very thin films on solid surfaces 
and for the measurement of optical constants of 
reflecting surfaces. The lower limit of film thick- 
nesses that can be studied by ellipsometry is at 
least an order of magnitude smaller than can be 
studied by other means such as interferometry. 
Artifacts such as those caused by vacuum in the 
case of electron microscopy are not encountered. 
Neither interferometry nor electron microscopy are 
adaptable, as is ellipsometry, to the study of films 
under liquids, and only ellipsometry will give the 
index of refraction of films of unknown thickness. 

The technique of ellipsometry is concerned with 
the measurement of changes in the state of polariza- 
tion of light upon reflection from a surface. For 
a clean reflecting surface the optical constants of 
the surface and the reflection coefficients of the 
system may be calculated from these changes. A 
thin transparent film on the reflecting surface causes 
additional changes from which the thickness and 
refractive index of the film may be determined. 

The principle of the ellipsometer has been de- 
scribed by several authors [e.g., 1-7] 2 and a bibli- 
ography of the theoretical contributions of many 
authors, starting with the original equations of 
Drude [1] is given by Winterbottom [2]. This paper 
will describe certain measurement techniques and 
an application of the exact solution of the algebraic 
equations of Drude to the measurement of optical 
constants of surfaces and thickness and refractive 
index of thin films covering the surfaces. A gen- 



i The work reported here was supported in part by the Army Research Office 
I Durham) and in part by the Bureau of Naval Weapons. 
2 Figures in brackets indicate the literature references at the end of this paper. 



orally applicable computational method is described. 
The complex and lengthy calculations required to 
solve the equations involved have been programmed 
for an electronic computer. Measurements carried 
out on chrome and gold surfaces with and without 
films will also be described. 

In addition, the Poincare sphere representation of 
the state of polarization of light will be developed 
since this is the most useful representation for the 
consideration of the effects occurring on reflection. 

2. Representations of Elliptically Polarized 
Light 

The state of polarization of elliptically polarized 
light may be described in many ways. For this 
purpose, consider a light wave traveling along the 
Z axis of a coordinate system. The electric vector 
of the wave is then given by: 



E x =a l cos (r+Si) 
E Y = d2 cos (r+<5 2 ) 



(i) 



where r=ult ]?co and v are the angular frequency 

and linear velocity of the light, respectively, and the 
total amplitude is the vector sum of a x and a 2 . The 
polarization can be described by the amplitudes, 
<Zi and a 2 , and the phase difference, d = 8 2 — 8 u of the 
components. As is well known, in a stationary plane 
whose normal is parallel to the Z axis the locus of the 
end of the electric vector is the ellipse shown in 
figure 1. 

The ellipse, however, may also be described in 
relation to coordinates X f and F' along the axes of 
the ellipse. Thus, the inclination <p of these coordi- 
nates and the semiaxes a and b of the ellipse also 
describe the polarization of the light. 
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Figure 1. The locus of the electric vector for elliptically 
polarized light in a plane normal to the direction of propa- 
gation. 

The coordinate axes ZFand X'Y' are both convenient for the description of 
the ellipse. 



It will be convenient to introduce auxiliary angles 
a and x defined by 



, a 2 

tan a=— 



tan x=±- 



(2) 
(3) 



The numerical value of tan x represents the ratio of 
the minor to major axes of the ellipse and the sign of 
X distinguishes the two senses in which the ellipse 
may be described. The angles x and y are called the 
ellipticity and azimuth of the light, respectively, and 
the representation may be made either in terms of 
&i, a 2 , and 8; a, 8, and the total amplitude; or <p, x, 
and the total amplitude. 

Another representation of the state of polarization, 
and one which leads quite naturally to the Poincare 
sphere is by parameters which all have the same 
physical dimensions, called Stokes parameters [8]: 



S 2 =2aia 2 cos 8 
S z =2a 1 a 2 sin 8. 



(4) 



The parameter S is proportional to the intensity of 
the wave and is related to the other parameters by 
the identity Sl = Sl+8%+Sl, so that the Stokes 
parameters are not all independent. 

It has been shown [9, 10] that the Stokes parameters 
are also given by 



Si=S cos 2x cos 2<p 
S 2 =S cos 2x sin 2p 
S 3 =S sm2 x . 



(5) 



Therefore, the state of polarization may be repre- 
sented by a point on a sphere of radius S by the 
spherical coordinates S , 2<p, 2%. This sphere is 
called the Poincare sphere. Alternatively, the 
Stokes parameters S ly S 2 , and S 3 may be used as 
Cartesian coordinates to describe the polarization. 
This representation is illustrated in figure 2. 

Since the intensity of the light is of secondary 
importance for use of the ellipsometer, the projection 
of the point on a sphere of unit radius is convenient. 
The coordinates are then given by 



*l = 



s 2 = 



-1 

1 — tan 2 a 

1 + tan 2 a 

2 tan a cos 8 
1+tan 2 a 

2 tan a sin 8 



1+tan 2 a 



cos 2% cos 2<p 
=cos 2x sin 2<p 
sin 2%. 



(6) 



Light of a given state of polarization is therefore 
represented by a point on the sphere, with the polar 
angles 2 <p and 2 x representing, respectively, twice 
the azimuth and twice the ellipticity of the light. 
For example, for plane polarized light, the ellipticity 
X is zero, and s 3 =0. Therefore, plane polarized light 
is represented by points on the equator of the 
Poincare sphere. Also, the ellipticity of circularly 
polarized light is 45°, whence Si = s 2 =0, 2x=±tt/2; 
thus, circularly polarized light is represented by the 
poles of the sphere. 

The Poincare sphere may be used to represent 
elliptically polarized light with respect to any 
physical axes, say X" and Y" at an angle (p" with 
the XY axis. The ellipticity x is not dependent on 
the choice of axes and the azimuth of the light is 
changed by <p". Thus the Si and S 2 axes must be 
rotated an angle 2 y" around the S 3 axis. 




>s 2 



Figure 2. 



Representation of polarized light on the Poincare 
sphere. 

The cartesian 



The point P represents light with ellipticity x and azimuth ty. 
coordinates are Stokes parameters (after Born and Wolf [9]) . 
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Figure 3. The intersection of the Poincare sphere and a plane 
normal to the Si axis at si = (l — tan 2 a)/(l-\-tan 2 a) for the con- 
sideration of the effect of a doubly refracting plate with re- 
tardation on light with ellipticity represented by P. 
The plate rotates P to L. The fast axis of the plate is along Si. 



The Poincare sphere is convenient for considera- 
tion of the effects of doubly refracting plates and 
reflection on the state of polarization of a light beam. 
For example, let the polarized wave represented in 
figure 1 pass through a doubly refracting plate of 
relative phase retardation 0. Choose the X and Y 
axes in figure 1 along the slow and fast axes of the 
plate. Since the azimuth of the fast axis of the 
plate is zero, it is represented by the positive S x 
axis. In this case there will be no change in the 
amplitudes a x and a 2 of the components of the light, 
and therefore no change in a or in s { as shown by 
eq (6), the only effect being to change the relative 
phase of the components by 0. Thus, upon passing 
through the plate the point representing the polariza- 
tion will remain on the curve representing constant 
«!, or on the intersection of the sphere and the plane 
s 1= constant. The intersection of the Poincare 
sphere by a plane perpendicular to the S x axis at s iy 

shown in figure 3, is a circle of radius — - — — ^— • 

l + tan 2 a 

The light incident on the plate is represented by 

the point P, and after passing through the plate is 

represented by the point L. The effect of the doubly 

refracting plate is seen to turn points on the sphere 

about the S x axis by an angle 0. If the plate had been 

oriented with its fast axis at an azimuth <p' with 

respect to the XY axes, the rotation by would have 

been about the line from the center of the sphere 

and the point 2 <p' on the equator. 

Reflection at a metal surface may be represented 

similarly. For reflection, the incident light wave is 

resolved into components in the plane of incidence 

and normal to the plane of incidence (the plane of 

the surface). The process of reflection introduces a 

phase difference A between these two components, 

and changes the ratio of their amplitudes by a 

factor tan \f/; that is, tan ^ is a measure of the rela- 



tive absorption of the two components. Thus, the 
ratio of the reflection coefficient for light polarized 
in the plane of incidence to that for light polarized 
in the plane of the surface is given (3) by 



p=— =tan \l/e j * 



(7) 



where p is the ratio of the reflection coefficients r p 
and r s , and \p and A are functions of the optical 
constants of the surface, the wavelength of the light 
used, the angle of incidence, and, for a film covered 
surface, the thickness and refractive index of the 
film. (See also eqs 40 to 46.) 

For consideration of the reflection process on the 
Poincare sphere, represent the light with respect to 
Cartesian coordinates with the X and Y axes in and 
normal to the plane of incidence and the Z axis in 
the direction of propagation. Then tan \p=aja v . 
The Stokes parameters of the light before reflection 
are given by eq (6) and after reflection are 



H 



1 — tan 2 a cot 2 \j/ 

1 + tan 2 a cot 2 ^ 

2 tan a. cot \p cos (5+ A) 

1+tan 2 a cot 2 yp 

2 tan a cot \p sin (5+A) 
1 + tan 2 a cot 2 \j/ 



(8) 



The reflection is represented on a Poincare sphere 
in figure 4. The point I v represents the intersection 





Figure 4. Representation of metallic reflection on the Ponicare 
sphere. 

The azimuth of the plane of incidence is represented by I P and the plane of the 
surface by L. In part a the view is along the Si axis. The point P is rotated t o 
L around the S\ axis by the retardation A of the surface. Eelative absorption 
of the components in and normal to the plane of incidence causes translation of 
L to L' along the great circle I p —L. The view in b is normal to this great circle . 
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of the axis Si with the sphere and has the Cartesian 
coordinate Si = l, s 2 = 0, s 3 =0. By eq (6) this cor- 
responds to zero ellipticity, X, and azimuth, <p, so it 
represents the orientation of the X axis, or the plane 
of incidence. Likewise, the point I s has Cartesian 
coordinates (—1,0,0), so it corresponds to zero 
ellipticity and an azimuth 7r/2. Therefore, the 
point I s represents the orientation of the Y axis or 
plane of the surface. The point P representing the 
polarization of the incident light is rotated around 
the Si axis through an angle A, the phase change 
produced by reflection, to the point L. The Stokes 
parameters of L are 



*i= 



1 — tan 2 a 



s 2 = 



S3-- 



l+tan 2 a 

2 tan a cos (8+ A) 
1+tan 2 a 

2 tan a sin (<5+ A) 
1+tan 2 a 



(9) 



In addition, the relative absorption of the com- 
ponents translates the point L to the point 27, the 
Stokes parameters for which are given by eqs (8). 

From eqs (8) and (9) we have s^/s f 2 = sys' 2 = tan (<5+ A). 
Points L and L f will, therefore, lie on the great circle 
given by the locus of points on the sphere with 
sjs 2 = tan (5+ A). This locus is shown in figures 4a 
and 4b and is given by the intersection of the sphere 
with a plane passing through the point L and normal 
to the axis Si. Thus the point L must be moved 
along the great circle defined above to point L f to 
represent the polarization of the reflected light. The 
position of L' on this great circle will now be derived. 
The Cartesian coordinates of I p and of I s are (1,0, 0) 
and (— 1, 0, 0), respectively, and the coordinates of L 
are given by eqs (9). Therefore, the length of the 
chords I P L and I S L are 



T j _ 2 tan a 

lp (1 + tan 2 a) 1 / 2 



I,L- 



(1 + tan 2 a) 1/2 



so 



IpL 
LL 



= tan a. 



(10) 



(11) 



The angle I P LI S is a right angle since it is inscribed 
in a semicircle. Therefore, the angle I V I S L is equal 
to a, and it is clear by eq (2) that the tangent of this 
angle is equal to the ratio of the amplitudes of the 
components of the incident light. 

The position of L' is determined by the angle 
I P I S L' . The tangent of this angle again gives the 
ratio of the amplitudes of the components of the 
reflected light, which is now equal to tan a/tan \p. 
Therefore, 



IM 



IsL 



7 =tan (I P I S L') = tan a/ tan \p 



(12) 



and the chord I V L' may be calculated as 

2 



IpL ' (l+cot 2 atan 2 ./<) 1/2 



Moreover, by eqs (11) and (12), we have 



tan 1^= 



T P L/I S L 
"l P L'/I s L' 



(13) 



(14) 



and, in general, the ratio of the chords between a 
point representing the state of polarization of a given 
light wave and any two diametrically opposed points 
on the sphere is equal to the ratio of the amplitudes 
of the components when the electric vector of the 
light is resolved along axes represented by these two 
opposed points. 

The Poincare sphere representation of the reflec- 
tion of light from a surface with known complex 
reflection coefficient tan \f/e jA may now be sum- 
marized. The point P, representing the state of 
polarization of the incident light is moved by the 
reflection process through an angle A on the sphere 
in a plane normal to the S x axis to the point L. The 
chords IpL and I S L are measured and tan a calculated 
by eq (11). Finally the angle I p I s L f is calculated 
from eq (12), or the chord I P L' is calculated from 
eq (13), thus determining the state of polarization 
of the reflected light. We shall show later how this 
representation ma}^ be used conveniently for the cal- 
culation of the effects occurring on reflection. First, 
however, we give a brief description of the instru- 
ment. 

2.1. Instrument 

Figure 5 shows the various components of an 
ellipsometer. Collimated monochromatic light, usu- 
ally the mercury green line (5460. 73 A), is used. 
The polarizer, a Glan-Thompson or a Nicol prism 
mounted in a graduated circle, serves to polarize the 
light emitted by the source. The compensator, also 
mounted in a graduated circle, is a birefringent pla'e 
usually of quarter-wave thickness; it is used to con- 



Light „ _ ... 

Source Colltmoting Lenses 



Collimated ~~"^£ 

Polychromatic ^^7^ 
Unpolarized "^ 



Monochromatic 
Light 




Elliptically 
Polarized 
Light 



Figure 5. The component parts of an ellipsometer. 
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vert the linearly polarized light into elliptically 
polarized light. The light incident upon tin 1 sample 
has an azimuthal angle and an ellipticity predictable 
from the settings of the polarizer and compensator. 
In general, the light will have its ellipticity and azi- 
muth changed by reflection from the sample. 

The "aperature" is an adjustable opening allowing 
a variation in the area of surface examined and the 
amount of light reaching the phototube. The 
analyzer, a second prism in a graduated circle, 
can be rotated until a minimum intensity is achieved, 
as indicated by the photometer. If some ellipticity 
exists in the light reaching the analyzer, extinction 
of the light cannot be obtained by rotation of only 
the analyzer. The polarizer and analyzer are then 
adjusted alternately to remove this ellipticity and 
obtain extinction. When this occurs, the light re- 
flected from the sample is plane polarized and can 
be thereby extinguished by the analyzer. The 
phototube and photometer permit determination of 
the null point with a high degree of sensitivity. 
For more accurate null point determinations graph- 
ical plots of the intensity of the transmitted beam 
versus angle of polarizer and analyzer may be used 

in]. 

2.2. Aline me nt 

Alinement of the ellipsometer is not too critical 
for the measurement of thickness and refractive 
index of thin films. Here the important quantity 
is the change in readings of the polarizer and ana- 
lyzer as compared to the values for the hare substrate. 
However, for obtaining accurate values of the optical 
constants of surfaces, alinement is* critical. More- 
over, during alinement certain confusing phenomena 
may occur, and it is worthwhile here to develop the 
theory of the alinement process. 

The procedure for alining the ellipsometer is typi- 
cal of that generally used in optical spectrometers 
except for the adjustment of the analyzer (A) and 
polarizer (P) scales. When alined, the scales for the 
polarizer and analyzer read zero when the planes of 
transmission of the prisms are parallel to the plane 
of incidence. This is accomplished by first adjusting 
the polarizer and analyzer prisms in their scales so 
tli at these scales differ in setting by 90° when the 
prisms are crossed. These scales, imagined as a unit, 
must then be rotated so that when the P and A scale 
read and 90° respectively, the planes of transmission 
of the two prisms are respectively in the plane of 
incidence and normal to the plane of incidence. That 
is, the coordinate system defined by the and 90° 
settings of both P and A must be made to correspond 
to the coordinate system defined by the plane of 
incidence and the plane of the surface. 

In principle, the adjustment is relatively simple. 
With the compensator removed and the polarizer 
and analyzer arms in the 'straight-through' position, 
the polarizer and analyzer prisms are adjusted 
in their respective scales until extinction is achieved 
with the scale readings differing by 90°. The arms 
arc then set for reflection from a metal surface. A 



minimum in the photometer reading is sought at 
which the scale readings differ by 90°. The planes 
of transmission of the two prisms should then be in 
the plane of incidence and normal to it. The prisms 
may now be rotated in their holders until the scale 
readings are ±180°, and 90 ± 180°, and the aline- 
ment is complete. 

In practice, however, certain subtle effects occur. 
After having first crossed the prisms when reflecting 
from the surface, a minimum may be achieved either 
by fixing P and adjusting A, or vice versa. If 
the light issuing from the polarizer is a ecu lately 
plane polarized, then the minimum achieved by 
either of these two means, and with the P and A 
scales crossed, occurs at a single value of P and A 
independent of the procedure followed. However, 
if the light issuing from the polarizer has some 
slight ellipticity due to some imperfection in the 
polarizer, then there are two positions at which the 
prisms are crossed and minimum photometer readings 
are obtained, depending upon whether the minimum 
is achieved by fixing P and adjusting A or vice 
versa. Moreover, the deepest minimum or most 
complete extinction of the reflected light now 
occurs at a third point, at which the scales are not 
crossed. Fortunately, with certain assumptions, 
even in this case alinement may be achieved quite 
accurately. We shall first give a theoretical explana- 
tion of the phenomena observed, and then give a 
procedure for alining the ellipsometer. 

We shall assume thai the polarizer produces light 
of a constant ellipticity x whose azimuth changes as 
the polarizer is rotated. The analyzer is taken to 
be perfect. It is assumed that the P and A prisms 
have been adjusted without reflection so that when 
minimum transmission is achieved the scale readings 
differ by 90°. The compensator is removed from 
the system. The process will be depicted on the 
Poincare sphere. We consider only ellipticity and 
azimuth much smaller than unity. In this case the 
Stokes parameters given by eq (6) are approximated 
by 

Sl =l 

s 2 = 2a cos<5 = 2(£> (15) 

s 3 = 2a sin<5 = 2x- 

Hence the state of polarization may now be repre- 
sented by the plane polar coordinates 2a and 8, 
and this corresponds to representing as a plane 
surface the portion of the sphere on the equator 
around the azimuth representing the plane of 
incidence, i.e., around I p . 

Under the above assumptions, the states of polar- 
ization of the light issuing from the polarizer all lie 
on the line 



2a sin 8 = 2x 



P= 



tan 5 



(16) 



(17) 



Any point P? with the polar coordinates 2a, 8, on this 
line is rotated and translated by the process of re- 
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flection to the point P' with coordinates 2a ', <5', with 
the relations 



2a' = 



2a 



tan \p 
8' = 5+A. 



(18) 



This process is illustrated in figure 6, where I—F 
represents the locus of points representing the inci- 
dent light and R—R' the reflected light. 

The state of polarization of the light after reflec- 
tion, represented by the point P', may be resolved 
into components along the plane of transmission of 
the analyzer and normal to it. The latter direction 
is represented in the figure by the point 2L. Under 
these conditions of very small azimuth and ellip- 
ticity, and by eqs (11) and (2), the amplitude of the 
light transmitted by the analyzer is proportional to 
the distance D between P' and 2L. 

By the triangle P', 2L, I p , the intensity of the 
transmitted light is proportional to 



D 2 =4L 2 - 



4x 2 



tan 2 \p sin 2 8 



SLx / cos A 
tan \p \tan 8 



-sin A 



(19) 



where 2L is to be considered an algebraic quantity, 
negative in figure 6. 

The process of taking readings may be accom- 
plished in two ways. First, the azimuth of P may 
be set and A adjusted for minimum transmission. 
This is equivalent to fixing 8 and adjusting L, i.e., 



d(£> 2 



dL 



=0, 



whence 



or 



L tan yp — X (cos A cos 8 — sin A) 
(^+|) tan y//=P cos A— x sin A (20) 



which is the equation of the locus of points obtained 
in this manner. 

Equation (20) is plotted in figure 7. The polarizer 
and analyzer will be crossed for only one point on 

7T 

this curve. For this point, P=A-\--> and 



P=- 



X sin A 



(21) 



cos A— tan \p 

and this crossed position is in general not at zero 

azimuth ( 8=- )> but is determined by the constants 

A and \p of the surface. 

However, the scales may also be adjusted by first 
setting A and moving P. In this case, the condition 
for a minimum in transmitted light is 



d8 



=0, L tan \f/ tan 5 cos A=x 



R 
I \ 


s 3 


\ 2X 1 

\ 2 L Y\ 


yi/8 
(11 


\^ y 


Ip 2P s 2 

^-2CX 
TAN f 

\ R 



Figure 6. The Poincare representation of metallic reflection 
for light of very small azimuth and ellipticity caused by an 
imperfect polarizer. 

The point along I— I' such as Pi is translated to a point on R—R' such as P' 
by the process of reflection. This is passed through as analyzer with plane of 
transmission at an azimuth X+7r/2. The amplitude of the light issuing from 
the analyzer is proportional to D. 



ABSOLUTE 
^^^ \ MINIMUM 


P / 
/ L = A+f 




XT /O 



Figure 7. The locus of points P and A representing extinction 
of reflected light for an imperfect polarizer and perfect 
analyzer. 

The upper line is achieved by fixing A and adjusting P for each setting, and 
the lower by fixing P and adjusting A for each setting. Note that the upper line 
passes through zero azimuth. 



and 



(-4+|) tan ^ cos A=P 



(22) 



which is different from eq (20). The point for which 
P and A are crossed is given bv 



P(l— tan t^cos A) = 0. 



(23) 



Since the quantity in parentheses is not in general 
zero, 

P=0. (24) 
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Thus it seems clear that the crossed position of the 
prisms achieved by setting A and adjusting P 
for a minimum will always be the true azimuth under 
these circumstances, i.e., an imperfect polarizer, but 
perfect analyzer. This is represented in figure 7, 
where we show both the line achieved by fixing P 
and adjusting A and the line achieved by fixing A 
and adjusting P. These curves were calculated from 
eqs (4) and (6) for A=135 deg and tan ^=K, and P 
and A are in the units of X. It will be noticed that the 
curve obtained by first fixing A passes through the 
zero azimuth. 

For comparison, a curve obtained experimentally 
is shown in figure 8, where the numbers indicate 
meter readings and hence are proportional to the 
transmission. Since the axes of figure 8 are inverted 
with respect to figure 7, the elliptic! ty of the polarizer 
must be negative. However, the similarity is strik- 
ing. The value of x computed from this figure and 
eq (21) is -0.18 deg. 

It will be clear from a consideration of eqs (20), 
(22), and (12), that it is best to use a polarization 
azimuth near the plane of incidence rather than near 
the plane of the surface, since tan \p is generally less 
than unity, and this has the effect of magnifying the 
angular difference between P and A±w/2 by an 
amount 1/tan \p. If the polarizer were set near the 
plane of the surface the angular spread between P 
and A~\r -k\2 would be decreased by an amount tan \p, 
thus decreasing experimental precision. 

The condition for the deepest minimum or best 
extinction is 

d(D 2 )_d(D 2 ) 



dL 



d<5 







and is thus given bv the simultaneous solution of 
eqs (20) and (22). this gives 



— x 




(25) 



ABSOLUTE MINIMUM 



Figure 8. The experimental counterpart of figure 7. 

The P and A coordinates of the intersection of the line P=A—ir/2 and the 
"fixing 1" line represent I lie azimuth of the plane of incidence and the plane of 
the surface, respectively. 



P= 



tan A 



By eqs (17) and (26) 



or 



tan A=— tan 8 

5+A = 7T 



(26) 



(27) 



so that the deepest minimum occurs when the re- 
flected light is linearly polarized, as would be ex- 
pected. This minimum is shown in figures 7 and 8 
by the intersection of the two lines. The absolute 
minimum determined by alternate adjustment of P 
and A is the same as this point of intersection within 
experimental error. It will be noticed that this point 

is not on the line P=A+— However, if x = 0, then 

the deepest minimum occurs at zero azimuth for 
both P and L, and thus either method of alinement 
may be used if the light is perfectly plane polarized. 
Thus, there appears to be only one unambiguous 
way of alining the ellipsometer when the light from 
the polarizer has ellipticity, namely, to seek that 

7T 

point at which both P=A±^ and minimum trans- 
mission is obtained when A is set and P adjusted. 

2.3. Alinement Procedure 

A step-by-step procedure for alinement of the 
ellipsometer is as follows: 

First, remove the quarter-wave plate. Lower the 
arms to the "straight through" position. Set the P 
scale to read zero or any convenient value. Adjust 
A until minimum transmission is achieved. The A 
scale will now in general read (P±90)+e. Rotate 
the A prism in its holder until e is zero. The P and 
A scales are now crossed when the prisms are crossed. 
The scales should track within ±0.02 deg throughout 
one revolution. 

With a metal surface set for reflection, raise the 
arms to an angle near the principal angle. Set A 
so that the plane of transmission of the analyzer is 
approximately in the plane of the surface. Adjust 
P until minimum transmission is achieved, and note 
the readings of the P and A scales. Move A by 
0.1 deg and again adjust P and note the values. If 
the meter reading is higher than for the previous set 
of values, adjust A by 0.1 deg in the other direction. 
If it is lower, continue in the same direction. Take 
a series of readings in this manner, changing A by 
0.1 deg and adjusting P for minimum transmission 
at each setting of A, making sure that the settings 
encompass the lowest meter reading. 

Plot the values of P and A and determine the 
readings of the P and A scales at which P=A±90°. 
At this point the plane of transmission of the P 
prism is in the plane of incidence and the plane 
of transmission of the A prism is in the plane of the 
surface. However, the scale readings will in general 
not be some multiple of 90°. Let the values of P 
and A at this point be, for example, 0° + e' and 90°+ e 
respectively. The quantity e' is now the amount by 
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which both scales are displaced from the true zero 
azimuth. Set A at the value obtained from the 
graph. Adjust the P prism in its holder until e is 
zero by turning the prism a small amount in its 
holder and again adjusting the P scale for minimum 
transmission. This will require a series of adjust- 
ments. (It is imperative that the A scale and prism 
not be moved at this stage and that the adjustment 
be done in this way. If, for example, the P scale 
were held and A adjusted, the ellipsometer would be 

alined to the value of P=A±- obtained by fixing P 

and adjusting A, which has been demonstrated to be 
incorrect.) Having adjusted e' to be zero, remove 
the reflecting surface and lower the arms to the 
"straight through" position. Set the P scale at 
zero and adjust the A scale for minimum. The A 
scale will now read 90+ e . Turn the A prism in its 
holder until e' is zero. The ellipsometer P and A 
scales are now alined. 

This procedure eliminates alinement errors due to 
ellipticity produced by the polarizer. The effect of 
ellipticity in the analyzer on alinement has not been 
determined but is expected to be small since the 
analyzer is used for extinction of the light. Thus, 
any ellipticity by the analyzer after extinction of the 
light will have no effect on the photometer reading. 

2.4. Determination of A and \p From Ellipsometer 
Scale Readings 

By eq (7) the reflection of light from a surface is 
characterized by the complex reflection coefficient 
tan \pe 3A , and hence by the two quantities \p and A. 
It will later be shown how the refractive index of a 
surface and the thickness and index of films on a 
surface are calculated from values of A and \f/. 
However, the determination of A and \p from P and 
A readings merits some discussion, for this is not 
always a perfectly direct matter. 

For any given surface there is a multiplicity of 
polarizer, analyzer, and compensator scale settings 
that produce extinction by the analyzer of the 
reflected light from the surface, and it becomes 
somewhat of a problem to determine the values of 
A and \p from these various readings. In order to 
explain how these numerous readings arise and how 
A and \f/ may be computed from them it is well to 
keep two facts in mind: (a) all azimuthal angles are 
measured positive counter-clockwise from the plane 
of incidence when looking into the light beam, and 
(b) the compensator, which may be set at any azi- 
muth, is generally set so that its fast axis is in an 
azimuth of ±t/4. The present discussion is for an 
instrument such as that shown in figure 5 with the 
compensator before reflection. If the compensator 
is placed after reflection, suitable corrections will 
have to be made. 

The various readings fall into four sets called 
zones, two with the fast axis of the compensator set 
at 7r/4, numbered 2 and 4, and two with it set at 
— t/4, numbered 1 and 3. In each zone there is 
one independent set of polarizer and analyzer read- 
ings, making four independent sets of P and A 



readings in all. However, since both analyzer and 
polarizer may be rotated by t without affecting the 
results, there are 16 polarizer and analyzer settings 
falling into four independent zones. Since the 
compensator may also be rotated by t without 
affecting the results, there are 32 possible sets of 
readings on the ellipsometer. 

Rather than calculating A and \p directly from the 
P and A values, it is useful to calculate three other 
quantities, p, a p , and a s from the P and A values, 
p being related to the P readings, a p related to the 
A readings in zones 1 and 4, and a s related to the 
A readings in zones 2 and 3. For a perfect quarter- 
wave plate these are related to A and \p by the eq (2) 



\}/=a 



a,. 



(28) 
(29) 



If the compensator is not a perfect quarter-wave 
plate, the relationships are 



— tan A=sin 8 cot 2p 
tan 2 i£=tan a p tan a s 



(30) 
(31) 



where 8 is the relative retardation of the com- 
pensator. If the retardation 8 is near 7r/2, then, to 
first order, eqs (28) and (30) are the same. However, 
it is now necessary to distinguish between a p and a s . 
The relationship of p and a p and a s to the P and A 
readings observed on the ellipsometer is best derived 
by a consideration of the process as represented on 
the Poincare sphere, and since this has been 
adequately treated by Winterbottom [2], we shall 
merely list the results. The meanings of p, a p , 
and a s in the four zones are as follows: 

Zone 1. The fast axis of compensator is at — 7r/4. 
The polarizer plane of transmission makes an angle 
of -\-p with the plane of incidence. The analyzer 
plane of transmission makes an angle of -\-a p with 
the plane of incidence. 

Zone 2. The fast axis of the compensator set 
at +7r/4. The polarizer plane of transmission 
makes an angle of — p with the surface (an angle 
of ir/2—p with the plane of incidence). The 
analyzer plane of transmission makes an angle of 
-\-a s with the plane of incidence. 

Zone 3. The fast axis of compensator is at — 7r/4. 
The polarizer plane of transmission makes an angle 
of -\-p with the surface (an angle of ^+7r/2 with 
the plane of incidence). The analyzer plane of 
transmission makes an angle of — a s with the plane 
of incidence. 

Zone 4-- The fast axis of compensator is at +7r/4. 
The polarizer plane of transmission makes an angle 
of — p with the plane of incidence. The analyzer 
plane of transmission makes an angle of — a p with 
the plane of incidence. 

The meanings of p, a p , and a s are now clear. The 
first is the angle between the polarizer plane of trans- 
mission and either the plane of incidence or the 
plane of the surface, while a p and a s represent the 
angle between the analyzer plane of transmission 
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and the plane of incidence, being labeled either a v 
or a s according to whether p is the angle between 
the polarizer plane of transmission and the plane 
of incidence or the plane of the surface, respectively. 
With these relationships, it is a simple matter to 
compute the relationship of P and A to p and a p 
or a s , and hence to A and \p, by means of eqs (28) 
and (29), or (30) and (31). These relationships are 
summarized in table 1. 

Table 1. Relation of P and A readings to p, a p , and a s 



Zone 


Compensator 


P 


A 


1 


-tt/4 


P 

p+ir 

P 

P + TT 


a p 

Op 

a p +7r 

flp + TT 


3 


-tt/4 


P+f/2 
P+3tt/2 
P+ir/2 
p+3ir/2 


7T— a s 
ir— a.< 
27r— a« 
2ir— a« 


2 


+W4 


r/2-p 

37T/2-P 

ir/2 - p 
Zir/2-p 


a g 

tta-j- IT 


4 


+»/4 


n-p 
2-rr-p 

JT-P 
2ir — p 


IT— dp 
IT — dp 

2-rr — dp 
2k -dp 



For the purposes of this table, it has been assumed 
that the P, A, and compensator scales have been 
adjusted to read angles with respect to the plane of 
incidence as previously defined. If the instrument 
is not arranged in this manner, it is necessary to 
adjust the readings before using table 1. 

The table gives the 16 possihle readings for the 
two settings of the compensator. When working 
with a completely unknown surface it is still not a 
simple matter to identify the P and A readings with 
one of the entries on this table and hence to com- 
pute j), a p , or a s . As a first step in this identification 
it is worthwhile to consider the ranges for A ii\\d\p 
and hence for p, a p , or a s . For all surfaces, including 
multiple reflections, 



In fact, in the large majority of cases 

o<*<? 



(32) 



(33) 



and from eq (29) the range of a p and a s is the same. 
For the special case of reflection from a dielectric 
surface at an angle of incidence smaller than the 
Brewster angle, A=—ir. For all other cases 



0<A<2tt 
which gives as a range for p, 



TT . 37T 

-i<p<x 



(34) 



(35) 



When working with a completely unknown surface, 
and particularly if multiple reflections are used, it is 
best to take a complete set of 16 readings for the 
two settings of the compensator. Then, with the 
relations (33) and (35), and the further condition that 
the values of p, a p , and a s calculated from these 
readings must be approximately the same, it is a 
simple matter to identify each of the scale readings 
with one of the possibilities outlined in the table, 
To continue to do this when the surface is known 
would be redundant and one reading in each zone is 
sufficient. We have found that the average of one 
reading in each of two zones with the same compen- 
sator setting gives approximately the same result as 
the average of one reading in each of all four zones. 

2.5. Differences Among the Zones 

In practice, the values of p, a p , and a, found from 
the P and A readings in the various zones are not 
identical. A typical set of values is shown in table 2. 
Differences of as much as 3 deg occur in the values of 
p and somewhat less in the values of a p and a s . On 
the other hand, it will be noticed that the average of 
the observed values of p in zones 1 and 3 is the same 
as the average of the value obtained in zones 2 and 4. 
This strongly implies the following equations 

Pi=p + 8 p i =p + 5 f 
p z =p—5 p2=p — b' 



(30) 



where 5 and <5' are error terms. These error terms are 
not constant from run to run, but are usually 1 to 
L.5 deg as in table 2. (In this particular set, 5 and 
8' are approximately the same, This is not true in 
general). They are not caused by misalinement or 
any other readily detectable cause, and their source 
remains obscure, However, among very many 
experiments, not a single one has been found where 
the averages from zones 1 and 3 did not check the 
average 4 from zones 2 and 4 within experimental error. 

Table 2. Values of P and A for a typical chrome slide in air 
in all four zones 



Zone 


P 


A 


P 


dp 


a s 


1 
2 
3 

4 


19.71 
L67. 17 
113. 00 
160. 20 


31. 04 
31.73 

147. 75 
148. 52 


19.71 

22. S3 

23. 00 
19. 80 


31.04 
31.48 


31.73 
32. 25 


Average of four zones 
Average of zones 1 & « 




21. 34 
21. 36 
21. 32 


31.63 


L 


31. 65 


L _ 


31. til 













The situation for a p and a s is not quite so clear-cut. 
If the compensator is not a perfect quarter-wave 
plate, then differences in a v and a s are to be expected. 
In fact, it may be shown that for a compensator with 
retardation 6 near t/2, 



-a s =(6—T/2) cos 2p sin 2^ 



(37) 



where, for the purposes of this calculation, we have 
used \p=a. Again, an inspection of the table indi- 
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cates that the average of a p and a s from zones 1 and 3 
is equal to the average of a p and a s from zones 2 and 
4 which again strongly implies 



(Li=a p -\-8 
a 3 =a s —5 



a 2 =a s —5' 



(38) 



where the a's with the numerical subscript are ob- 
served values in the various zones. Many experi- 
ments indicate that these equations seem to be 
followed quite faithfully. The values of 5 and 5' in 
eqs (38) are, in general, different from those in 
eqs (36). 

It will be observed that the equations are not 
independent, and they cannot be solved for a p and 
a s directly from the observed values of a. Fortu- 
nately, for a small value of a p —a s , and for \[/ not 
too close to zero it may be shown that 



(tan a p tan a 5 ) 1/2 = tan 



(a P +a s ) 



(39) 



to first order in a p —a s . Thus rather than using eq 
(31) for computing yp, the above equation is used. 
This requires a compensator not too different from 
quarter-wave, and independent measurements indi- 
cate that the compensator used in this work has a 
retardation which differs from w/2 by about 1 deg. 

Because the averages from zones 1 and 3 check so 
closely the averages from zones 2 and 4 for both p 
and a, measurements are usually taken in zones 1 
and 3 only and averaged. 

2.6. Surfaces 

The substrates most easy to work with are those 
with high reflectance, such as metals. Both smooth- 
ness and flatness are factors that must be considered 
in the selection of a substrate. Irregularities that 
are small compared to the dimensions of the light 
beam, which is commonly of the order of 1 mm 2 , are 
averaged and do not affect the results [4]. Long- 
range regularity (flatness) is desirable if more than 
one location on a specimen is to be studied. An indi- 
cation of the regularity is obtained by determination 
of p and a at several points on the specimen. Meas- 
urements on 1-in. long steel gage blocks of 0.09 /x in. 
tolerance varied only about 0. 1 ° in the polarizer readings 
from top to bottom, and less in analyzer readings. 
Acceptable slides prepared by shearing small rec- 
tangles from a large ferrotype plate varied about 0.2° 
for similar lengths and these make convenient and 
readily available surfaces for study. 

Various techniques have been used for the prepara- 
tion and cleaning of the substrate surface. The 
chrome ferrotype plate slides used in these investiga- 
tions were washed in warm distilled organic solvents 
to remove organic contamination. Since the result- 
ing surface was hydrophobic, the slides were further 
cleaned with warm chromic acid cleaning solution 
followed by several washings in warm distilled water. 
This treatment should not appreciably affect the 
character of the chromium-chromium oxide surface, 
but it did remove organic contamination, the surface 



now being hydrophilic. However, in less than 1 hr 
the surface became sufficiently contaminated to 
become hydrophobic, even when enclosed in a covered 
container. These precleaned slides were passed 
through a flame immediately before use to remove 
this contamination; this flaming of the slides restored 
the hydrophilic character to the surface. For studies 
under liquids, the slides were placed under the liquid 
while still warm from the flame. The cleaning of 
surfaces by flaming has been described by Patrick 
[12] and Bartell and Betts [13]. 

2.7. Cells 

Cells have been used for ellipsometer measure- 
ments in vacuum or gaseous environments [2, 14, 15] 
and under liquids [2]. The cell shown in figure 9 has 
been used in our studies for measurements under 
liquids. The specimen is situated on the base of the 
cell which is filled with liquid of known refractive 
index. The light beam enters and leaves the cell 
through optically flat windows. These windows are 
inclined at the angle of incidence $ with respect to the 
base of the cell so that the light passes through them 
at normal incidence. Therefore, reflection of light 
at the surface of the cell windows will be independent 
of its direction of polarization and the polarization of 
the light will not be changed. There should be no 
stresses in the glass sufficient to cause detectable 
birefringence; moreover, the inner and outer sides of 
each window should be parallel. They may be sealed 
to the cell body either by fusion or by epoxy resin; 
the epoxy seal has been used more successfully for 
our cells. The differences in polarizer and analyzer 
readings in air for a metallic reflecting surface out- 
side and inside the cell were no larger than 0.1° for 
our cells. 

Difficulties with photometer fluctuations have oc- 
curred as a result of convection currents caused by 
evaporation of liquid. A small area of liquid in con- 
tact with the air, a tightly fitting cover, and solu- 
tions filtered through fritted glass disks minimize 
this effect. 

The angle of incidence should generally be chosen 
to give the maximum sensitivity for the measure- 
ment of film thickness. For this purpose, sensitivity 
may be defined as the change of P reading or A 
reading with film thickness, t, i.e., dP/dt and dA/dt. 
One is then faced with the problem of selecting an 
angle of incidence such that these two quantities 
are a maximum. No general rules can be laid out 
for this selection, and the choice of angle of incidence 
will depend upon the particular substrate, film and 
surrounding medium. However, for the common case 
of an organic film on a chromium surface in an 
organic medium, the sensitivity for both P and A 
has been calculated for various angles of incidence 
and film thicknesses up to 1000 A. The results are 
shown in figures 10 and 11. 

The maximum sensitivity for A occurs at angles 
of incidence between 75 and 80 deg, and decreases 
with film thickness. For P, the maximum sensitiv- 
ity (the absolute value of dP/dt) occurs at an angle 
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Figure 9. A cell used for studies of surfaces under liquids. 
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Figure 10. The sensitivity, dA/dt, o/ A /or measuring thick- 
ness t, plotted as a function of angle of incidence. 

The parameters are those for a typical organic film on chrome; n x = 1.4268, 
712= 1.4700, m =3.14800-4.142797. 
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Figure 11. The sensitivity, dP/dt of V for measuring thickness t, 
for the same conditions as figure 10. 
Note the very low sensitivity at a film thickness of 700 A. 



of incidence of 70°, and decreases with increasing 
film thickness up to a thickness of about 700 A. 
At this thickness, the sensitivity in P is very small, 
making the thickness determination entirely depend- 
ent on A, which is relatively insensitive in this region. 
Thus, under these conditions, the accurate determi- 
nation of film thickness is difficult for films 600 to 
800 A thick. For films thicker than 800 a, the sen- 
sitivity in P again becomes usable, but the maximum 
occurs at angles less than 70 deg. 

For any film thickness, therefore, the best angle 
of incidence is always a compromise between the 
most sensitive angle for P and that for A. The 
selection will, in part, depend on the relative impor- 
tance of P or A to the determination of the thickness. 

These sensitivities apply only to the given specific 
conditions but the sensitivities for other conditions 
with a nonabsorbing film are expected to have similar 
behavior although specific values will be different. 

2.8. Methods of Computation 

A typical system for study by the ellipsometer 
consists of a film of index n 2 and thickness d on a 
reflecting substrate of index n :i immersed in a 
medium of index n\, as shown in figure 12. Let all 
media be isotropic and n x represent a real index 
of refraction, while n 2 and 7i 3 may be complex. 

Consider light incident at the boundary between 
the immersion medium and film. The cosine of 
the refraction angle is 



cos 



^=[ l -(S sin ^)] ,/2 - (40) 



The parallel and normal reflection coefficients for 
light incident at this boundary are: 



and 



p n 2 cos cpi—rii cos <p 2 

12 7l 2 COS <pi+^i COS <p 2 



_ n x cos <p\—n 2 cos <p 2 
~i cos <pi+n 2 cos <P2 



n* 



(41) 
(42) 



mmersion 
Medium 




Figure 12. Reflection from a film-covered surface. 
The substrate is a reflecting metal surface. 
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respectively. The reflection coefficients, rf 3 and r 23 
at the boundary between the film and substrate are 
given by similar expressions. 

The total reflection coefficients, R p and R s that 
include the contributions of reflections from lower 
boundaries are given [2] by: 



R p 



rf 3 exp D 



and 



1+^2^3 exp D 



Rs __ rU+rk exp D 
l+r s 12 r s 23 exp D 



(43) 



(44) 



where cos <p 3 values needed for these reflection 
coefficients are given by an expression similar to 
eq (40) and D represents the quantity 



D=—4cTrjn 2 cos <p 2 d 2 /X 



(45) 



where X is the w avele ngth of the light used, in 
vacuum, and ./ = V — 1. The ratio of the parallel 
and normal total reflection coefficients is defined 
as p: 

P=R*/R\ (46) 

This may be expressed in terms of the relative 
attenuation and phase shift of the parallel component 
with respect to the perpendicular component that 
occurs, represented by the azimuthal angle A, and 
relative phase shift \p by: 



p=tan yp exp (jA) 



(47) 



nr- 



as in eq (8). Thus p is determined from ellipsom- 
eter readings. 

The value of the complex index of a reflecting 
surface can be calculated from the equation 

[\ 4psiir>ri 1/2 
=fh tan ^[l- (p+1) rj (48) 

where ^ is the angle of incidence and p is determined 
from ellipsometry measurements on the base sub- 
strate. 

Several methods of determining the thicknesses of 
films on reflecting substrates from ellipsometry meas- 
urements are available [2-5]. When the indexes of 
the substrate and film are known, tables or graphs of 
A and \p may be computed from given values of d. 
This is accomplished by calculating values of p from 
eqs (40) to (46) and values of A and \p from eq (47). 
Figure 13 shows curves computed for three different 
values of film index n 2 . The numbers along each 
curve are thickness values of the film, d 2 . Experi- 
mental values of A and \p are then interpolated in 
such tables or related graphs to determine unknown 
thicknesses. 

However, it is usually more efficient to solve the 
equations directly from the thickness of a film. Sub- 
stituting eqs > (43), (44), and (46) in eq (47) and 
rearranging gives a quadratic of the form: 



CI (exp D) 2 +C 2 (exp D) + C 3 = 



(49) 
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Figure 13. A plot of A against \p for films of refractive index 
as shown, in a medium with a refractive index of 1.359. 

The numbers along the curves are film thicknesses in Angstrom units. From 
such a figure and an experimental point \p e , A«, both the refractive index and 
thickness of the film may be obtained (see text) . 

where C\, C 2 , and (7 3 are complex functions of the 
refractive indexes, angles of incidence, A and \p. 
For a given value of the coefficients eq (49) gives two 
solutions for exp D and a value of d may be calculated 
from each. Since the coefficients are complex, the 
film thicknesses calculated from this equation would 
also be expected to be complex. However, the cor- 
rect film thickness, d, must be a real number as it 
represents a real quantity. Therefore, the solution 
of the quadratic that yields a real film thickness is 
the correct solution. In practice, various experi- 
mental errors will result in both solutions yielding 
complex values for d. The thickness with the small- 
est imaginary component is selected as the correct 
solution; the real portion is taken as d and the imagi- 
nary part, dj, is taken as a relative measure of error. 

The real portion of d is then used to compute A and 
\p by eqs (43) to (47). As the imaginary component 
of d has been dropped, these values will differ from 
the experimental angles by amounts <5A and dxj/, and 
dj, 3 A, and 5^ are all measures of the experimental 
error. However, <5A and 5\f/ must be within the limits 
of experimental error of yp and A for the results to be 
valid. This is a more direct determination of the 
validity of an experiment than the magnitude of dj. 

If both the thickness and index of refraction of the 
film, n 2 , are not known, the equations cannot be 
solved for d and n 2 in closed form. For this case, a 
series of refractive indexes are assumed and a thick- 
ness is calculated from the experimental measure- 
ments. These calculations will result in error terms, 
described above, of different magnitudes. The 
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range of refractive indexes and related thicknesses 
within the experimental errors <5A and b\j/ are then 
chosen. The magnitude of the possible ranges of 
refractive indexes and thicknesses will depend upon 
the magnitude of the error and on several parameters 
such as rii, n 2 , n 3 , etc. 

Some of the principles described above are 
illustrated in figure 13, which shows theoretical A and 
\p values for a range of given film indexes and thick- 
nesses. An experimentally determined point (A e , 
yp e ) for a film of unknown n 2 falls near one of these 
curves. For an assumed film index of 1.5, a thick- 
ness will be calculated corresponding to the point 
(A c , \p c ). The error terms, 8A and 5^ are also illus- 
trated. The error terms are usually much smaller 
than indicated on this figure. 

The calculations discussed above can be com- 
puted manually, but the equations are complicated 
and the computation of any appreciable amount of 
data would he impractical. Several "Fortran" 
programs have been written in this laboratory for 
IBM 704 and 7090 computers to enable the com- 
putation of data from ellipsometer measurements. 
The programs are general enough to permit their use 
in most situations encountered. The following is a 
brief description of problems that can be solved with 
these programs. 

The refractive index of the substrate, n, 3 , may be 
computed directly from experimental readings on the 
bare substrate. Immersion media of different re- 
fractive indexes may be used for this determination. 

When a p and a s are known, the correct retardation 
of the compensator can be calculated. This retarda- 
tion may then be used in subsequent calculations. 
Thus, the ellipsometer data may be used even if the 
compensator is not a quarter-wave plate. 

Curves and tables of A \p, and reflection coefficients 
can be computed for a series of given film thicknesses 
and refractive indexes, n 2 , as shown in figure 13. 

For experimental values of A and \p, a thickness, <L 
and error terms, 5^ and 8 A, are computed for a given 
n 2 . If n 2 is unknown, a series of n 2 can be assumed 
and corresponding film thichness and error terms 
computed. The selection of the correct n 2 and 
thickness has been discussed and is illustrated in 
the next section. All the above calculations have 
been extended for multiple films and for multiple 
reflections. 

3. Applications 

Inasmuch as surfaces in air will usually have an 
adsorbed film of water or other contaminants, in 
order to determine the index of a substrate it is 
necessary either to measure the surface in a vacuum 
or in a medium with a refractive index identical to 
that of the adsorbed film (n { =n 2 ). 

In the work reported here chrome slides with and 
without a vacuum deposit of gold were flamed to 
remove adsorbed gases and immediately immersed 
in a liquid. Table 3 gives the refractive indexes of 
slides measured in air and measured under water. 
The differences in the refractive index of the metals 
calculated from measurements made in air and under 



Table 3. Measured refractive index, n = N — jK, of chrome 
and gold a in air and in water 



Immersion medium 


Chrome 


Gold 




N 


K 


N 


K 


Air 


2.954 

:<._>3i 

3. 005 
3.218 

2.950 
3.285 


4. 224 
4.354 

4. 274 
4. 302 

4. 258 
4.426 


0. 440 
. 528 

. 427 
. 540 




2.134 


Water 


2. 105 


Air 


2. No 


Water 


2.100 


Air 




Water. 













a Slides had been flamed prior to examination. 

water are probably caused by the presence of an 
adsorbed film of water on the slides when measured 
in air. Such a film would not be observed, of course, 
when the measurements were made under water, 
assuming- the refractive index of the film to be the 
same as that of water. It is presumed, therefore, 
that the measurements made under water yielded 
the correct indexes for these slides. 

The measurements described above were also 
repeated for immersion media of different refractive 
indexes. The optical constants of the metal calcu- 
lated under these conditions are «;iven in table 4. 



Table 4. Refractive indexes, N — jK, calculated for chrome. 
slides immersed in various liquids immediately after flaming 



<p 


Immersion medium 


n 


Retractive index ■< 




N 


K 


Degret 
60 


methanol.. . ... 


L.329 

1.337 
1.337 
1.301) 
L.360 
L.426 
1.498 

1.00 
1.00 


3.212 

3.410 
3. 278 
3. 275 
:>. 296 
3. 254 
3.344 

•2. 899 
2.918 


4.335 


60 


water ... . . 


4.358 


7(1 


water . .. 


4.371 


60 


acetone 


4.363 


70 


acetone . .. ._ . 


4.398 


70 


cyclohexane ._ 


4.396 


70 


toluene.— .. 


4.437 


60 


air b _ _ 


1. 170 


7ii 


air'- 


4.198 









■•<■ All refractive Indexes are averages of two or more sets of measurements, ex- 
cept for methanol. 

i> Slides measured in air had been exposed to air from 1 to 24 hr. 



The slides measured under liquids were flamed and 
i in mediately immersed in the liquid before a film 
could form on them from the air. These slides are 
therefore not expected to have any adsorbed film on 
them, so that their true refractive indexes are 
measured. If some film did remain on them, its 
refractive index would probably be near that of the 
immersion liquid, causing only a small error in the 
measured refractive index of the slide. Hence the 
rerfactive indexes of the slides under the various 
liquids are all approximately the same, and inde- 
pendent of the liquid and the angle of incidence. 
The variation is due mainly to variations among the 
individual slides used for the different measurements. 
The refractive index measured in air shows a 
difference due to neglect of the films that is larger 
than the variations among the individual slides. 
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This film is expected to be adsorbed water and, 
perhaps, other gases and is expected to have an 
index near that of water. By first measuring slides 
in air and then measuring them under water, it is 
therefore possible to calculate a thickness of the 
adsorbed layer in air, assuming the refractive index 
of the adsorbed film to be the same as that of liquid 
water. Such results are shown in table 5, where 
<5A and 5^ are the differences obtained as described 
above. The reproducibility of the results is seen to 
be quite good, although the actual value of the 
thickness may be somewhat in error due to the 
assumption of the refractive index value. 

Table 5. Thickness of adsorbed layer a on slides in air 





Slide 


<P 


A 


5A 


<A 


b^ 


d 




number 


















Deqree 


Deqree 


Deqree 


Deqree 


Deqree 


A 


Chrome 


76 


50 


163. 44 


0.00 


39.51 


0.04 


23 




76 


60 


152. 56 


.02 


3(3. 33 


. 11 


24 




76 


70 


132. 52 


.06 


31.85 


.13 


26 




75 


70 


132. 72 


.08 


31.57 


. 16 


27 




74 


70 


131. 78 


.04 


31.83 


.08 


24 




79 


70 


132. 04 


.06 


31.76 


. 12 


23 




76 


80 


85.04 


.06 


28.15 


.18 


26 




25 






Gold 


72 


70 


90.36 


0. 06 


40. 53 


0. 15 


2 




73 


70 


91.96 


.00 


40. 98 


.02 


3 




74 


70 


88.38 


. 12 


39.89 


. 26 


5 




77 


70 


88.36 


.06 


40.39 


.16 


4 




4 







1 Calculated as water 
T=22 ° C, 
RH= 50 percent. 



The determination of both the thickness and 
refractive index of a barium fluoride film vacuum 
evaporated onto a chrome surface is given as an 
example of the measurements and calculations 
described earlier. The index, n^ of the substrate 
was first determined from the reading taken on the 
bare chrome surface, then the film was vacuum 
deposited on the surface and the angles A and \f/ 
were measured for the film-covered surface. A 
series of refractive indexes are assumed for the film, 
and thickness and the corresponding error terms, 
<5A and 5^, are calculated for each index by the 
methods given earlier. The results are given in 
table 6. The reproducibility of the readings is 
±0.1° for A and ±0.05° for \p and the maximum 
experimental error is expected to be about ±0.2° 
for A and ±0.1° for \j/. It will be observed from 
table 6 that the index 1.456 gives zero error terms 
8A and 5\p, and is, therefore, the best fit with the 
experimental reading. However, n 2 values from 
1.454 to 1.458 also give error terms within the limits 
of experimental error, and are the range of possible 
indexes of the film. The corresponding thicknesses 
are 654 and 665 A. 

The actual range of possible thicknesses is greater 
than this range as is illustrated in figure 14. This 
figure is an enlargement of the center portion of 
figure 13. Curves for additional refractive indexes 



Table 6. Thickness 



of evaporated 
chrome slide a 



barium fluoride film on 

b 



A 


8A 


•A 


5^ 


n 2 


d 


Degree 


Degree 


Degree 


Degree 




A 


133. 96 


3.64 


36.30 


0.60 


1.420 


753 




2.64 




.40 


1.430 


726 




1.58 




.22 


1.440 


701 




1.08 




.14 


1.445 


688 




0.58 




.07 


1.450 


676 




.28 




.03 


1.453 


668 




.18 




.02 


1.454 


665 




.10 




.01 


1.455 


662 




.00 




.00 


1.456 


660 




.10 




.01 


1.457 


657 




.20 




.02 


1.458 


654 




.28 




.03 


1.459 


651 




.38 




.04 


1.460 


647 




.70 




.38 


1.475 


453 



a Immersedin acetone, ?ii = 1.359, at an angle of incidence of 60°. The refractive 
index of the slide was m =3.316—4.383;'. 
b Film prepared by Frank E. Jones, National Bureau of Standards. 



have been added. The dashed lines represent the 
limits of experimental error for the point (A e , \j/ e ). 
The true values of n 2 and d must fall within these 
limits. A range of possible thicknesses for each 
refractive index rather than only a single selected 
value of thickness for each possible refractive index 
is obtained. The range of possible thicknesses is 
thus approximately 630 to 680 A. 

These methods do not determine the index of the 
film as accurately when the film is very thin. This 
can be observed in figure 13 by the close proximity 
of the curves for small thicknesses. Since the 
experimental error is independent of film thickness, 
the range of refractive indexes and thicknesses that 
would be included in the rectangle of experimental 
error described above is increased. When the 
refractive index is known, the thickness can be 
accurately determined even for very thin films, 
except that as the film index approaches that of the 
medium, the sensitivity in determining the film 
thickness decreases. This can be seen in figure 13 
by the decreasing sensitivity to thickness of A and \p 
as the refractive index of the film, n 2} approaches the 
index of the medium, 1.359. 




FiCxURE 14. An enlargement of figure 13, showing how error 
in Ae andxpe leads to an uncertainty in the calculated refractive 
index and thickness of a film. 
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4. Summary 

The use of the ellipsometer for the measurement 
of the thickness and index of refraction of very thin 
films is described and illustrated by examples. The 
representation of elliptically polarized light by means 
of the Poincare sphere is developed in detail since it 
is necessary for a complete understanding of the 
changes in the polarization of light in the ellipsometer. 

Some unusual effects in alinement have been 
observed and explained as due to the light from the 
polarizer having a slight ellipticity. A method of 
alinement that is not affected by the ellipticity is 
given. 

Chrome ferrotype plates have been found to be a 
convenient substrate for study of thin films. Means 
of preparing and cleaning the plates are discussed. 
Also, cells used to study films under liquids are 
described. 

In order to study a film by ellipsometry, the 
reflection coefficient of a bare substrate is first 
measured and the complex refractive index of the 
substrate computed from the reflection coefficient. 
A film is then deposited on the substrate and the 
reflection coefficient of the combination is measured. 
If the index of the film is known, the thickness of 
the film may be computed. In one method, the 
reflection coefficient of the substrate with a film is 
computed for many thicknesses Oi the film. Then 
the measured reflection coefficient is interpolated in 
the calculated results to determine the thickness 
of the film. However, in this paper the required 
equations are solved to give the thickness of the 
film directly in terms of the measured reflection 
coefficients so that interpolation is not required. 
These measurements and calculations are applied 
to determining the thickness of an adsorbed film on 
ferrotype plates. The index of refraction of this 
film was assumed to be that of water. 

Both the index of refraction and thickness of a 
film may be calculated from the complex reflection 
coefficient of the film-substrate combination. The 
procedure is to assume a series of refractive indexes 
and compute a corresponding thickness for each 



index. If the refractive index chosen is not the 
exact index of the film, and if the measurements are 
not sufficiently precise, the calculated thickness of 
the film is complex. The imaginary part of the 
computed complex thickness is taken as a relative 
measure of error either of the assumed refractive 
index or of the original measurement. The film 
refractive index with its corresponding thickness 
that yields the smallest error terms is taken as the 
best fit. This calculation is illustrated for a barium 
fluoride film. 
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